\(\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx\) [223]
Optimal result
Integrand size = 46, antiderivative size = 37 \[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=-\frac {B \cos (e+f x) (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n}{f}
\]
[Out]
-B*cos(f*x+e)*(a-a*sin(f*x+e))^m*(c+c*sin(f*x+e))^n/f
Rubi [A] (verified)
Time = 0.08 (sec) , antiderivative size = 37, normalized size of antiderivative = 1.00,
number of steps used = 1, number of rules used = 1, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.022, Rules used = {3049}
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=-\frac {B \cos (e+f x) (a-a \sin (e+f x))^m (c \sin (e+f x)+c)^n}{f}
\]
[In]
Int[(a - a*Sin[e + f*x])^m*(c + c*Sin[e + f*x])^n*(B*(m - n) + B*(1 + m + n)*Sin[e + f*x]),x]
[Out]
-((B*Cos[e + f*x]*(a - a*Sin[e + f*x])^m*(c + c*Sin[e + f*x])^n)/f)
Rule 3049
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_) + (d_.)*sin[(e_.
) + (f_.)*(x_)])^(n_.), x_Symbol] :> Simp[(-B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/(f*
(m + n + 1))), x] /; FreeQ[{a, b, c, d, e, f, A, B, m, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && EqQ
[A*b*(m + n + 1) + a*B*(m - n), 0] && NeQ[m, -2^(-1)]
Rubi steps \begin{align*}
\text {integral}& = -\frac {B \cos (e+f x) (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n}{f} \\
\end{align*}
Mathematica [A] (verified)
Time = 1.31 (sec) , antiderivative size = 37, normalized size of antiderivative = 1.00
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=-\frac {B \cos (e+f x) (c (1+\sin (e+f x)))^n (a-a \sin (e+f x))^m}{f}
\]
[In]
Integrate[(a - a*Sin[e + f*x])^m*(c + c*Sin[e + f*x])^n*(B*(m - n) + B*(1 + m + n)*Sin[e + f*x]),x]
[Out]
-((B*Cos[e + f*x]*(c*(1 + Sin[e + f*x]))^n*(a - a*Sin[e + f*x])^m)/f)
Maple [A] (verified)
Time = 9.80 (sec) , antiderivative size = 38, normalized size of antiderivative = 1.03
| | |
| method | result | size |
| | |
| parallelrisch |
\(-\frac {B \left (-a \left (\sin \left (f x +e \right )-1\right )\right )^{m} \left (c \left (1+\sin \left (f x +e \right )\right )\right )^{n} \cos \left (f x +e \right )}{f}\) |
\(38\) |
| | |
|
|
|
|
[In]
int((a-a*sin(f*x+e))^m*(c+c*sin(f*x+e))^n*(B*(m-n)+B*(1+m+n)*sin(f*x+e)),x,method=_RETURNVERBOSE)
[Out]
-B/f*(-a*(sin(f*x+e)-1))^m*(c*(1+sin(f*x+e)))^n*cos(f*x+e)
Fricas [A] (verification not implemented)
none
Time = 0.26 (sec) , antiderivative size = 37, normalized size of antiderivative = 1.00
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=-\frac {{\left (-a \sin \left (f x + e\right ) + a\right )}^{m} {\left (c \sin \left (f x + e\right ) + c\right )}^{n} B \cos \left (f x + e\right )}{f}
\]
[In]
integrate((a-a*sin(f*x+e))^m*(c+c*sin(f*x+e))^n*(B*(m-n)+B*(1+m+n)*sin(f*x+e)),x, algorithm="fricas")
[Out]
-(-a*sin(f*x + e) + a)^m*(c*sin(f*x + e) + c)^n*B*cos(f*x + e)/f
Sympy [F]
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=B \left (\int m \left (- a \sin {\left (e + f x \right )} + a\right )^{m} \left (c \sin {\left (e + f x \right )} + c\right )^{n}\, dx + \int \left (- n \left (- a \sin {\left (e + f x \right )} + a\right )^{m} \left (c \sin {\left (e + f x \right )} + c\right )^{n}\right )\, dx + \int \left (- a \sin {\left (e + f x \right )} + a\right )^{m} \left (c \sin {\left (e + f x \right )} + c\right )^{n} \sin {\left (e + f x \right )}\, dx + \int m \left (- a \sin {\left (e + f x \right )} + a\right )^{m} \left (c \sin {\left (e + f x \right )} + c\right )^{n} \sin {\left (e + f x \right )}\, dx + \int n \left (- a \sin {\left (e + f x \right )} + a\right )^{m} \left (c \sin {\left (e + f x \right )} + c\right )^{n} \sin {\left (e + f x \right )}\, dx\right )
\]
[In]
integrate((a-a*sin(f*x+e))**m*(c+c*sin(f*x+e))**n*(B*(m-n)+B*(1+m+n)*sin(f*x+e)),x)
[Out]
B*(Integral(m*(-a*sin(e + f*x) + a)**m*(c*sin(e + f*x) + c)**n, x) + Integral(-n*(-a*sin(e + f*x) + a)**m*(c*s
in(e + f*x) + c)**n, x) + Integral((-a*sin(e + f*x) + a)**m*(c*sin(e + f*x) + c)**n*sin(e + f*x), x) + Integra
l(m*(-a*sin(e + f*x) + a)**m*(c*sin(e + f*x) + c)**n*sin(e + f*x), x) + Integral(n*(-a*sin(e + f*x) + a)**m*(c
*sin(e + f*x) + c)**n*sin(e + f*x), x))
Maxima [F]
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=\int { {\left (B {\left (m + n + 1\right )} \sin \left (f x + e\right ) + B {\left (m - n\right )}\right )} {\left (-a \sin \left (f x + e\right ) + a\right )}^{m} {\left (c \sin \left (f x + e\right ) + c\right )}^{n} \,d x }
\]
[In]
integrate((a-a*sin(f*x+e))^m*(c+c*sin(f*x+e))^n*(B*(m-n)+B*(1+m+n)*sin(f*x+e)),x, algorithm="maxima")
[Out]
integrate((B*(m + n + 1)*sin(f*x + e) + B*(m - n))*(-a*sin(f*x + e) + a)^m*(c*sin(f*x + e) + c)^n, x)
Giac [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 8948 vs. \(2 (37) = 74\).
Time = 51.92 (sec) , antiderivative size = 8948, normalized size of antiderivative = 241.84
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=\text {Too large to display}
\]
[In]
integrate((a-a*sin(f*x+e))^m*(c+c*sin(f*x+e))^n*(B*(m-n)+B*(1+m+n)*sin(f*x+e)),x, algorithm="giac")
[Out]
-(B*cos(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)
+ 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 2*pi*m*floor(-1/8*sgn(4*t
an(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*
x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*t
an(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)
+ 1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x +
e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4*pi*n)*e^(-m*log(2) - n*log(2) +
m*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*ta
n(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 +
abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan
(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2
- 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e
) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1
/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2))*abs(a)/(tan(f*x + e)^2*tan(
1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*t
an(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 +
8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^
2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*
e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) +
4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(
f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x +
1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(
1/2*f*x + 1/2*e)^2 + 1)))*tan(-1/4*pi*m*sgn(2*a*tan(1/2*f*x + 1/2*e)^4 - 4*a*tan(1/2*f*x + 1/2*e)^3 + 4*a*tan(
1/2*f*x + 1/2*e) - 2*a)*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)
) - 1/4*pi*n*sgn(2*c*tan(1/2*f*x + 1/2*e)^4 + 4*c*tan(1/2*f*x + 1/2*e)^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2*c)*sgn
(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)) + 1/2*pi*m*floor(f*x/pi +
e/pi + 1/2) + 1/2*pi*n*floor(f*x/pi + e/pi + 1/2) - 1/4*pi*m*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x
+ 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2
+ 4*c*tan(1/2*f*x + 1/2*e)))^2*tan(1/2*f*x + 1/2*e)^2 + 2*B*e^(-m*log(2) - n*log(2) + m*log(sqrt(2)*sqrt(abs(
4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2
*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan
(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*
tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(
1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2
*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2
+ 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2))*abs(a)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan
(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 +
8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) +
2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan
(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 +
abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*ta
n(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x +
1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*
x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1))
)*sin(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) +
4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 2*pi*m*floor(-1/8*sgn(4*tan
(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x
+ 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan
(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) +
1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)
^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4*pi*n)*tan(-1/4*pi*m*sgn(2*a*tan(1
/2*f*x + 1/2*e)^4 - 4*a*tan(1/2*f*x + 1/2*e)^3 + 4*a*tan(1/2*f*x + 1/2*e) - 2*a)*sgn(4*a*tan(1/2*f*x + 1/2*e)^
3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(2*c*tan(1/2*f*x + 1/2*e)^4 + 4*c*tan
(1/2*f*x + 1/2*e)^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2*c)*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e
)^2 + 4*c*tan(1/2*f*x + 1/2*e)) + 1/2*pi*m*floor(f*x/pi + e/pi + 1/2) + 1/2*pi*n*floor(f*x/pi + e/pi + 1/2) -
1/4*pi*m*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sg
n(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)))*tan(1/2*f*x + 1/2*e)^2
- B*cos(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)
+ 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 2*pi*m*floor(-1/8*sgn(4*t
an(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*
x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*t
an(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)
+ 1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x +
e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4*pi*n)*e^(-m*log(2) - n*log(2) +
m*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*ta
n(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 +
abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan
(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2
- 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e
) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1
/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2))*abs(a)/(tan(f*x + e)^2*tan(
1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*t
an(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 +
8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^
2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*
e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) +
4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(
f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x +
1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(
1/2*f*x + 1/2*e)^2 + 1)))*tan(-1/4*pi*m*sgn(2*a*tan(1/2*f*x + 1/2*e)^4 - 4*a*tan(1/2*f*x + 1/2*e)^3 + 4*a*tan(
1/2*f*x + 1/2*e) - 2*a)*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)
) - 1/4*pi*n*sgn(2*c*tan(1/2*f*x + 1/2*e)^4 + 4*c*tan(1/2*f*x + 1/2*e)^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2*c)*sgn
(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)) + 1/2*pi*m*floor(f*x/pi +
e/pi + 1/2) + 1/2*pi*n*floor(f*x/pi + e/pi + 1/2) - 1/4*pi*m*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x
+ 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2
+ 4*c*tan(1/2*f*x + 1/2*e)))^2 - B*cos(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(
f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) +
5/8) + 2*pi*m*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) +
4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e
)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)
^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan
(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4
*pi*n)*e^(-m*log(2) - n*log(2) + m*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x +
e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x
+ e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1
/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f
*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x +
1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 -
8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) +
2))*abs(a)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqr
t(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)
^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(
f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x +
1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f
*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan
(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*
tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x +
1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)))*tan(1/2*f*x + 1/2*e)^2 - 2*B*e^(-m*log(2) - n*log(2)
+ m*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*
tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2
+ abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*t
an(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)
^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2
*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x +
1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2))*abs(a)/(tan(f*x + e)^2*ta
n(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2
*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2
+ 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e
)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/
2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)
+ 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*ta
n(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x
+ 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + ta
n(1/2*f*x + 1/2*e)^2 + 1)))*sin(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)
^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 2
*pi*m*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f
*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e)^2*tan(
1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*t
an(1/2*f*x + 1/2*e) + 2) + 1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x
+ 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4*pi*n)*t
an(-1/4*pi*m*sgn(2*a*tan(1/2*f*x + 1/2*e)^4 - 4*a*tan(1/2*f*x + 1/2*e)^3 + 4*a*tan(1/2*f*x + 1/2*e) - 2*a)*sgn
(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(2*c*tan(1/
2*f*x + 1/2*e)^4 + 4*c*tan(1/2*f*x + 1/2*e)^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2*c)*sgn(4*c*tan(1/2*f*x + 1/2*e)^3
+ 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)) + 1/2*pi*m*floor(f*x/pi + e/pi + 1/2) + 1/2*pi*n*flo
or(f*x/pi + e/pi + 1/2) - 1/4*pi*m*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f
*x + 1/2*e)) - 1/4*pi*n*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)
)) + B*cos(2*pi*n*floor(-1/8*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*
e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 2*pi*m*floor(-1/8*sgn(
4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2
*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) + 5/8) + 1/4*pi*n*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 +
8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) +
2) + 1/4*pi*m*sgn(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x
+ e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2) - 1/4*pi*m - 1/4*pi*n)*e^(-m*log(2) - n*log(2
) + m*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4
*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2
+ abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*
tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e
)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/
2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(f*x + e)^2*tan(1/2*f*x
+ 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 - 8*tan(1/2*f*x + 1/2*e) + 2))*abs(a)/(tan(f*x + e)^2*t
an(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + tan(1/2*f*x + 1/2*e)^2 + 1)) + n*log(sqrt(2)*sqrt(abs(4*tan(f*x + e)^
2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2
+ 8*tan(1/2*f*x + 1/2*e) + 2)*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*
e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1
/2*e) + 2)*tan(f*x + e)^2 + abs(4*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e
) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2)*tan(1/2*f*x + 1/2*e)^2 + abs(4*t
an(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + 8*tan(f*x + e)^2*tan(1/2*f*x + 1/2*e) + 4*tan(f*x + e)^2 + 2*tan(1/2*f*
x + 1/2*e)^2 + 8*tan(1/2*f*x + 1/2*e) + 2))*abs(c)/(tan(f*x + e)^2*tan(1/2*f*x + 1/2*e)^2 + tan(f*x + e)^2 + t
an(1/2*f*x + 1/2*e)^2 + 1))))/(f*tan(-1/4*pi*m*sgn(2*a*tan(1/2*f*x + 1/2*e)^4 - 4*a*tan(1/2*f*x + 1/2*e)^3 + 4
*a*tan(1/2*f*x + 1/2*e) - 2*a)*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x +
1/2*e)) - 1/4*pi*n*sgn(2*c*tan(1/2*f*x + 1/2*e)^4 + 4*c*tan(1/2*f*x + 1/2*e)^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2
*c)*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)) + 1/2*pi*m*floor(f
*x/pi + e/pi + 1/2) + 1/2*pi*n*floor(f*x/pi + e/pi + 1/2) - 1/4*pi*m*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(
1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1
/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)))^2*tan(1/2*f*x + 1/2*e)^2 + f*tan(-1/4*pi*m*sgn(2*a*tan(1/2*f*x + 1/2*e)^4
- 4*a*tan(1/2*f*x + 1/2*e)^3 + 4*a*tan(1/2*f*x + 1/2*e) - 2*a)*sgn(4*a*tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f
*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(2*c*tan(1/2*f*x + 1/2*e)^4 + 4*c*tan(1/2*f*x + 1/2*e)
^3 - 4*c*tan(1/2*f*x + 1/2*e) - 2*c)*sgn(4*c*tan(1/2*f*x + 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2
*f*x + 1/2*e)) + 1/2*pi*m*floor(f*x/pi + e/pi + 1/2) + 1/2*pi*n*floor(f*x/pi + e/pi + 1/2) - 1/4*pi*m*sgn(4*a*
tan(1/2*f*x + 1/2*e)^3 - 8*a*tan(1/2*f*x + 1/2*e)^2 + 4*a*tan(1/2*f*x + 1/2*e)) - 1/4*pi*n*sgn(4*c*tan(1/2*f*x
+ 1/2*e)^3 + 8*c*tan(1/2*f*x + 1/2*e)^2 + 4*c*tan(1/2*f*x + 1/2*e)))^2 + f*tan(1/2*f*x + 1/2*e)^2 + f)
Mupad [B] (verification not implemented)
Time = 12.44 (sec) , antiderivative size = 37, normalized size of antiderivative = 1.00
\[
\int (a-a \sin (e+f x))^m (c+c \sin (e+f x))^n (B (m-n)+B (1+m+n) \sin (e+f x)) \, dx=-\frac {B\,\cos \left (e+f\,x\right )\,{\left (-a\,\left (\sin \left (e+f\,x\right )-1\right )\right )}^m\,{\left (c\,\left (\sin \left (e+f\,x\right )+1\right )\right )}^n}{f}
\]
[In]
int((B*(m - n) + B*sin(e + f*x)*(m + n + 1))*(a - a*sin(e + f*x))^m*(c + c*sin(e + f*x))^n,x)
[Out]
-(B*cos(e + f*x)*(-a*(sin(e + f*x) - 1))^m*(c*(sin(e + f*x) + 1))^n)/f